Transport in partially equilibrated inhomogeneous quantum wires 
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We study transport properties of weakly interacting one-dimensional electron systems including 
on an equal footing thermal equilibration due to three-particle collisions and the effects of large-scale 
inhomogeneities. We show that equilibration in an inhomogeneous quantum wire is characterized 
by the competition of interaction processes which reduce the electrons total momentum and such 
which change the number of right- and left-moving electrons. We find that the combined effect of 
interactions and inhomogeneities can dramatically increase the resistance of the wire. In addition, 
we find that the interactions strongly affect the thermoelectric properties of inhomogeneous wires 
and calculate their thermal conductance, thermopower, and Peltier coefficient. 

PACS numbers: 71.10.Pm, 73.23.-b 
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I. INTRODUCTION 

Transport properties of low-dimensional systems have 
been subject of intensive research work over the last two 
decades. One of the fundamental discoveries that has 
driven the field was the observation of conductance quan- 
tization in ballistic quantum wires and quantum point 
contactsji^ It was found that conductance exhibits a 
staircase- like dependence on the electron density with the 
universal step. The understanding of this phenomenon 
follows already from the single-electron picture, which 
predicts for the conductance of a one-dimensional single 
channel-clean wire,— 
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The physical origin of conductance plateaus at certain 
gate voltages was associated with a fixed number of occu- 
pied electronic subbands, each supplying one quantum of 
conductance 2e^ /h. Within the same approach of nonin- 
teracting particles both charge and energy are carried by 
electronic excitations. This results in the universal rela- 
tion between electric and thermal conductances, known 
as the Wiedemann-Franz law K = {n^/3e^)TG. The 
thermal conductance of noninteracting electrons is thus 
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In addition to G and K two thermoelectric coefficients 
of the electron gas are usually of great interest. These 
are thermopower S, which relates an induced voltage 
drop across the wire to applied temperature gradient, 
and Peltier coefficient 11 connecting electric and heat cur- 
rents. These two coefficients are connected by an Onsager 
relation 11 = ST. In the absence of interactions the ther- 
mopower and Peltier coefficients are exponentially small 



at low temperatures T ^ fi. (Here fj, is the chemical 
potential.) The reason for such strong suppression of 
thermopower and Peltier coefficients is the partial can- 
cellation between heat currents carried by particles with 
energies ^ + e and /x — e. Only the absence of electronic 
states below the bottom of the band prevents 11 and S 
from vanishing exactly. 

The remarkable success of the simple single-electron 
picture in describing the quantization of conductance 
and explaining the temperature dependence of thermo- 
electric coefficients is attributed to the fact that quan- 
tum wires are always connected to two-dimensional leads, 
where interactions between electrons do not play a sig- 
nificant role. Even though the interactions in the wire 
are usually not weak, i.e., e^/hvp ^ 1, where vp is the 
Fermi velocity, it has been shown within the so-called 
Luttinger-liquid model of one-dimensional electrons that 
the interactions inside the wire do not affect conductance 
quantization4r"— It is no surprise then that a number 
of recent experiments f^^— that revealed deviations from 
the perfect quantization, ((T|), in low-density wires, at- 
tracted a great deal of theoretical attentionJ^"— These 
deviations often take the form of a shoulder-like fea- 
ture, which develops at finite temperature just below the 
first quantized plateau, around 0.7 x 2e^//i. At present 
there is no consensus on the theoretical interpretation of 
this phenomenon. However, it is generally accepted that 
electron-electron interaction effects should be involved in 
explaining these experimental observations. 



In a number of recent publications 
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properties of one-dimensional conductors were reconsid- 
ered focusing on the physics which lies beyond an ideal 
Luttinger-liquid model. In particular, when studying the 
temperature dependence of the corresponding kinetic co- 
efficients Refs. [3ll - [39l emphasized one fundamental as- 
pect of interactions, namely the role of physical processes 
that lead to equilibration of electrons inside the wire. 
It should be emphasized that equilibration is absent in 
an ideal Luttinger liquid since bosonic elementary ex- 



citations of the latter have infinite lifetime, thus there 
is no relaxation towards equilibrium in these systems, 
no matter how strong the interactions are. In higher- 
dimensional systems equilibration at low temperatures is 
primarily provided by pair collisions of electrons. These, 
however, do not provide relaxation in one-dimensional 
systems. This is due to the conservation laws for momen- 
tum and energy which severely restrict the phase space 
available for scattering. As a result, pair collisions in 
ideal one-dimensional wires can occur with a zero mo- 
mentum change or an interchange of the two momenta, 
leaving the distribution function unaffected. The leading 
equilibration mechanism thus involves collisions of more 
than two particles. For a weakly interacting system, 
it is then natural to assume that equilibration is pro- 
vided by three-particle scattering processesJ,?::^! This, 
of course, also relies on the additional assumption that 
other degrees of freedom, which can absorb energy and 
momentum from electrons (phonons, for example) can 
be ignored. This assumption is acceptable in many cases 
since electron-phonon coupling constant is typically much 
smaller than that due to the electron-electron interac- 
tions. 

In practice, long one-dimensional structures are 
strongly prone to inhomogeneities inevitably present due 
to the nearby gates or charged dopants underlying the 
wire. However, most preceding works studied effect 
of equilibration on transport assuming uniform (clean) 
wires. The notable exceptions include Refs. '32 and [33 
where smooth inhomogeneities were accounted for while 
assuming full equilibration of the electronic system. The 
purpose of the present work is to study effects of inho- 
mogeneities on transport properties of partially equili- 
brated quantum wires. We focus our attention to the 
situation where the scale of inhomogeneities b is much 
larger than the electron Fermi wavelength, b ^ Xp, such 
that backscattering of electrons from the inhomogeneities 
is negligible. In this case only interactions (three-particle 
collisions) may interrupt direct flow of the electron liquid 
and convert (backscatter) some right-moving electrons 
into left-moving ones. 

Since non-uniform systems are no longer trans- 
lationally invariant, there is an additional scatter- 
ing mechanism, which can relax electron momen- 
tum without changing the number of right- and left- 
moving particles. '^^'■^■^ We find that equilibration due 
to three-particle collisions and inhomogeneity-induced 
momentum-nonconserving scattering compete with each 
other. An interplay between these two effects leads to a 
very interesting picture of the electronic transport and 
results in temperature-dependent corrections to the wire 
resistance and thermoelectric coefficients. 

The paper is organized as follows. In SeclUwe discuss 
the general structure of the electron distribution function 
and transport in a clean one-dimensional wire. We briefly 
mention the resulting transport coefficients for this case, 
which were recently reported in Ref. [S^. In Sec. IIIII we 
develop a general formalism which enables us to treat 



equilibration of the electron system due to three-particle 
collisions and inhomogeneity-induced scattering on equal 
footing. The central part of our work is Sec. IIVI where 
we apply this theory to the calculation of transport co- 
efficients in one-dimensional inhomogeneous wires. From 
our general expressions we recover the results known for 
the short and long uniform wires, and also consider sev- 
eral experimentally relevant simple models of inhomoge- 
neous case. We summarize our results in Sec.|Vl Supple- 
mentary appendices accompany some technical aspects 
of our calculations. 



II. TRANSPORT IN UNIFORM QUANTUM 
WIRES 



In the absence of interactions, left- and right-moving 
electrons inside the wire are at equilibrium with the reser- 
voir they originated from. If a voltage bias V and/or tem- 
perature difference AT is applied between the reservoirs, 
then corresponding equilibria differ from each other, giv- 
ing rise to a particular form of the nonequilibrium dis- 
tribution function inside the wire. This distribution de- 
pends on the direction of motion of electrons and for the 
right- and left-moving particles is controlled, respectively, 
by the left and right lead. 
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cp — p /2m is the energy of an electron with mo- 



Here €p 

mentum p and 9{p) is the unit step function. The dif- 
ference between the chemical potentials (temperatures) 
in the leads is equal to the voltage (temperature differ- 
ence) applied to the wire /i; — /x^ = eV {Ti — T^ = AT) . 
Using the distribution function (U) at AT = one can 
find electric current / ~ GV with the conductance of 
noninteracting electrons G = (2e^//i)(l -f e~''/"^)~^, thus 
recovering, Eq. ([T]), up to an exponentially small cor- 
rection. The same distribution, Eq. (|3]), provides ther- 
mal conductance, Eq. ([2]), and thermoelectric coefficients 
IV = TS= (/x/e)e-^/^, consistent with Eq. ©. 

In the presence of interactions ballistic propagation of 
electrons through the wire may be interrupted by col- 
lisions with other electrons. As a result of these colli- 
sions, some electrons change their direction of motion 
thus losing memory of the lead they originated from. 
Such backscattering processes modify the electron distri- 
bution function which is then no longer given by Eq. (j4|) . 
It is important to realize that the effect of electron col- 
lisions on the distribution function depends strongly on 
the length of the wire. Indeed, electrons traverse short 
wires relatively fast, such that interactions do not have 
time to change distribution, Eq. dU, considerably. On 
the other hand, in the limit of very long wire one should 
expect full equilibration of left- and right-moving elec- 
trons into a single distribution, even in the case of weak 
interactions. 



For the Galilean invariant system one can easily infer 
the electron distribution function in a fully equilibrated 
state. Indeed, viewed from a reference frame moving with 
the drift velocity Vd = I/ne (where I is the electric cur- 
rent and n is the electron density) the electron system is 
at rest and must be described by the equilibrium Fermi 
distribution. Performing a Galilean transformation back 
into the stationary frame of reference this distribution 
takes the form, 
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where the chemical potential /ieq and temperature Tcq 
inside the equilibrated wire are, in general, different from 
fJ-i(r) and T;(r). 

At zero temperature, T = Teq = 0, the distributions, 
Eqs. (|H) and (O, coincide, provided /i;(r) — Meq ± VdPF, 
where pp — ■Khn/2 is the Fermi momentum of the sys- 
tem. At non-zero temperature the distribution function, 
Eq. ([5]) , of electrons inside the equilibrated wire is slightly 
different from the distribution, Eq. (jj]), supplied by the 
leads. The mismatch between the two distribution func- 
tions results in additional resistance, reducing the con- 
ductance of noninteracting electrons to^ 
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This result is universal since it was obtained without 
making any specific assumptions regarding the process 
of equilibration. The corresponding derivation relied 
uniquely on the analysis of conservation laws for energy, 
momentum, and particle number. It is applicable as long 
as wire length L exceeds certain equilibration length l^q 
such that the distribution, Eq. ([S|), is already established. 
The exact definition of £eq is model specific and depends 
on the interaction between electrons. Quite generically, 
however, it can be argued that this length is exponen- 



tially large at low temperature 
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e^/^. The expo- 



nential scale can be understood from the mechanism of 
equilibration'^^ which is also discussed later in the text. 
The thermal conductance of fully equilibrated wire is 
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This result can be understood simply from the structure 
of the distribution function ^ . One should recall that 
thermal conductance is defined under the condition that 
electric current / = envd vanishes. Thus Vd = and 
the distribution, Eq. ([5]), takes the form of the standard 
Fermi-Dirac distribution. Due to its symmetry p — >■ —p 
the heat current carried by electrons vanishes, regardless 
of the temperature bias AT applied to the wire. One 
therefore finds that in an infinitely long wire the ther- 
mal conductance is zero. One should note, however, that 
thermal conductivity is finite, "^^ i.e., at L — > cx3 the ther- 
mal conductance scales a.s K oc 1/L. 



Unlike the electric and thermal conductances, ther- 
mopower and Peltier coefficients are significantly en- 
hanced by equilibration effects. Specifically, 11 grows 
from the exponentially small value, Eq. ([3]), for short 



wires L <^ i^q to 
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in fully equilibrated (long) wires, L ^ £oq- 

A more careful treatment of the equilibration effects is 
required in wires of intermediate length L ~ ^cq, where 
the electron distribution function is only partially equili- 
brated)^ As we already mentioned, in the case of weakly 
interacting electrons the leading mechanism of equilibra- 
tion is provided by three-particle collisions. At low tem- 
peratures one should consider two types of such colli- 
sions. The strongest scattering events involve three par- 
ticles near the Fermi level (for example, one left mover 
that scatters off two right movers such that all particles 
preserve their direction of motion). These collisions are 
relatively fast, and the corresponding scattering length £t 
scales as a power law of temperature. However, these col- 
lisions alone can not establish the distribution ([5|) since 
they conserve the number of right- and left- moving parti- 
cles. The other important three-particle collisions involve 
backscattcring of, say, a right-moving electron into a left- 
moving one. This backscattcring occurs near the bottom 
of the band and provides equilibration between the chem- 
ical potentials of right and left movers, thus establishing 
the distribution, Eq. ([5]). Since at low temperatures the 
probability to find an empty state at the band bottom is 
exponentially small, the corresponding relaxation process 
is very slow, and equilibration length is large, ^eq oc e^' "^. 
Let us consider now a segment of the wire, whose 
length Ai is small compared to the equilibration length 
^cq but large as compared to £t, namely, £t ^ Ai ^ ^eq- 
This condition implies that typical electron with energy 
near the Fermi level passes through the segment without 
backscattcring so that the distribution, Eq. ([S]), cannot 
be established. On the other hand, AL is already suf- 
ficiently large for electrons to experience other multiple 
collisions which allow momentum and energy exchange 
between right- and left-moving electrons. Under these 
conditions, the electron distribution function in the seg- 
ment achieves a state of partial equilibration, in which 
the numbers N^ and N^ of the right- and left-moving 
electrons are conserved independently. The form of this 
distribution can be obtained from the general statisti- 
cal mechanics argument by maximizing the entropy of 
electrons while preserving N^^^' , total energy and mo- 
mentum of the system f2§, 
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Here T is the effective temperature, parameter u has di- 
mension of velocity and accounts for the conservation 
of momentum in electron collisions, and [i^^^> are the 




FIG. 1: [Color online] Inhomogeneous quantum wire con- 
nected adiabatically to two-dimensional leads and biased by 
the small voltage fii— ^r = eV and/or temperature difference 
Ti — Tr = AT. We assume that the spatial scale b of the 
inhomogeneities is large compared to the Fermi wavelength, 
and thus electrons do not experience backscattering from in- 
homogeneities. In this case only three-particle equilibration 
processes may interrupt the flow of right-moving electrons 
and convert some right movers into left movers. Wavy lines 
represent heat exchange Q^ between right movers and left 
movers. 



chemical potentials of the left- and right-moving parti- 
cles. The distribution, Eq. ©, smoothly interpolates 
between the regimes of no equilibration, Eq. (j4]), and 
that of full equilibration, Eq. ([S]). In the absence of tem- 
perature difference, AT — 0, the unperturbed distribu- 
tion, Eq. (m, is obtained from Eq. ^ by setting u = 
and identifying the chemical potentials with those in the 
leads: /i^ — fii and ^^ = Hr- [Here and throughout the 
paper we use l{r) to denote left (right) lead while L{R) 
denote left (right) movers.] The fully equilibrated dis- 
tribution, Eq. (O, is obtained from Eq. (|9]) by setting 
A/x = jjL^ — M^ =0. In this case the electric current 
is expressed as / = enu, which identifies parameter u 
as the drift velocity Vd- We should emphasize here that 
since the distribution, Eq. (|9l), is applicable for the seg- 
ment of the wire outlined above then all four parameters 
T (x),m(x), and /i^'^(x) defining fp are, in principle, co- 
ordinate dependent. 

The implications of the distribution, Eq. (jH]), for 
the transport coefficients of partially equilibrated clean 
wires were discussed in Ref. [S^. In the following 
we generalize the above picture of electronic transport 
in one-dimensional wires accounting for possible non- 
uniformities of the system. 



III. TRANSPORT IN INHOMOGENEOUS 
QUANTUM WIRES 

A. Boltzmann equation 

Consider an inhomogeneous quantum wire of length 
L, connected by ideal reflectionless contacts to nonin- 
teracting leads and biased by a small voltage V and/or 
temperature difference AT, see Fig.[T] If the spatial vari- 
ations related to inhomogeneities occur on a length scale 
h much larger than the Fermi wavelength \p, electrons do 



not suffer any backscattering. Since the physical picture 
of equilibration in one-dimensional wire can be readily 
understood at the level of weakly interacting electrons 
we restrict our attention to this case and describe the 
system in the framework of kinetic equation. In this case 
electron distribution function f{t,x,p) obeys the Boltz- 
mann equation. 



dtf + vpd^f ~ d^U{x)dpf = I{f} , 



(10) 



where static potential U{x) accounts for inhomogeneities 
of the wire and T{/} conventionally stands for the col- 
lision integral. We are interested in the steady-state 
regime when distribution function f{t,x,p) does not de- 
pend explicitly on time, and thus set dtf = 0. It will be 
also convenient to split the distribution f{x,p) into two 
parts corresponding to the right and left movers, 

/(x,p) = e{p)f''ix, ep{x)) + 0(-p)/^(x, ep{x)) , (11) 

and express it as the function of energy ep(x) — p^/2m + 
U{x) for the given momentum p. Kinetic Eq. (|10p should 
be supplemented by the boundary conditions at the ends 
of the wire that are controlled by the leads. 
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where fii ^ ii + eV, j-ir — M ^-^d Ti = T + AT, 
Tr = T. [We use shorthand notation for the distribution 
function of right movers at the left lead f^'{l,ep{l)) = 
f^{x — Q,ep{x = 0)) and left movers at the right lead 
f^{r,ep{r)) — f^{x = L,ep{x = L)).] The parametriza- 
tion, Eq. ([TT|) , is especially useful since owing to the sim- 
ple algebraic relation. 
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the inhomogeneity-related term drops out from the left- 
hand side of the kinetic equation, except for the resid- 
ual contribution S{p)d^Uix)[f^{x,Uix)) - f^ix,U{x))] 
at p = 0. For noninteracting electrons the mismatch 
between distribution functions f-^/^{x,U{x)) of right 
and left movers is exponentially small at the bottom 
of the band. In addition, even this small discontinu- 
ity is smeared by inter-electron scattering responsible 
for equilibration.^^ It is thus safe to take f^{x,tp{x)) — 
f^{x, £p{x)) forp = and we get then instead of Eq. (|10l) . 

0{p)vpdj''{x, ep{x)) + 0(-p)«p9,/^(x, ep{x)) = T{/} . 

(14) 
As the first step of our general analysis, we demon- 
strate now with the help of kinetic Eq. p^ that devi- 
ations in electric and thermal conductances from their 
noninteracting values [Eqs. ([TJ and ([2])] are ultimately 
related to the rate of change in the number of say right- 
moving electrons N^ and heat exchange rate Q^ between 
right movers and left movers. 



B. Conservation laws 

The rate of change in the number of right movers N 
due to electron coUisions is obtained from the colhsion 
integral I{/} upon integration over positive momenta 
and wire length, 
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dx / dpX{f} , 



(15) 



where the coefficient 2 stands for two spin projections. 
Owing to the Boltzmann Eq. (1141) N^ can be equivalently 
presented in terms of the distribution function of right- 
moving electrons as 



iV« 
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We can integrate this expression by parts by noticing 
that 



U{x) 
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d_ 
dx 

and approximating in the following the distribution func- 
tion of right movers by unity at the bottom of the band, 
which is correct up to exponentially small terms at low 
temperatures f^{x, C/) « 1 - 0[e-'^^'-^^/'^]. This would 
give then 

nR = fir) ~ fil) + '^[U{r) - U{1)] , (18) 

where we used standard definition for the currents of 
right /left movers. 
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In Eq. ([T8)l the incoming current f{f) of right movers at 
X = is known since it is controlled by the distribution of 
noninteracting electrons in the left lead, Eq. ((T^ . How- 
ever, the outgoing j^{r) is not known because the dis- 
tribution function of right-movers varies along the wire 
as a result of scattering. It is convenient to exclude this 
unknown from Eq. (J18l) by noticing that the total current 
j in the wire does not depend on position due to the con- 
servation of the number of electrons and can be written 
as j = j^{r) + J^(?')- This results in 

N''=3- [f{r)+f{l)] + l[U{r) - U{1)] . (20) 

The benefit of performing this step is that now both 
currents j^{l) and j^{r) are controlled by the noninter- 
acting leads whose distribution functions are given by 
the boundary conditions [Eq. (fT2|)]. Furthermore, since 
N^ and j vanish in the absence of applied bias, we can 
exclude the [/-dependent contribution from Eq. (I^IH) by 



subtracting from it j^(r) + J^(Olv=o " 
This leads to 



;[C/(r)-C/(0]. 



iV«=j-[j^(0-J^(Ol^^o 



(21) 



The difference between currents of right movers at the 
left boundary with and without bias in Eq. (j2ip can be 
found with the help of distribution function (IT^ . Indeed, 
after a simple calculation 



valid up to corrections small as e '^'"^, we find 
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where / = ej. This result can be thought of as a 
generalization of Landauer formula for interacting one- 
dimensional systems. In the noninteracting limit N^ — 
and we recover G — I/V — 2e^//i while a finite N^ would 
lead to a change in the conductance. Equation (|22t was 
derived earlier for uniform (clean) wires<^i^ We have 
shown here that it remains intact even in the case of in- 
homogeneous wires. 

We now repeat the above calculation for the energy 
change of right-movers E^ induced by electron collisions. 
The latter is obtained from the collision integrallj/} by 
multiplying it by ep and then integrating over positive 
momenta and the wire length 
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(23) 



It can be equivalently rewritten in terms of f^{x,e) by 
making use of the Boltzmann Eq. (ITil) , 

^«= / dx deed^.fix,e). (24) 

" Jo Ju(x) 

After integration by parts, similar to Eq. ([T7|. one finds 

E'' = jE{r)~j§{l) + l[U'{r)-U\l)] , (25) 

where we used the usual definition for energy currents of 
right/left movers, 

Je i^) ^ -f^ dpO{±p)vpepf{x,ep). (26) 

Conservation of energy ensures that the energy current 
JE is constant along the wire. By using it at the right 
end JE = Jsi^) + il;(''): '^^ can exclude unknown j§{r) 
from Eq. (|25|) in analogy with Eq. ([20]) . In addition, 
since E^ = and je = without the bias, we subtract 
jhir)+Jmvlo' = UuHr) - U^l)] to exclude the 



^/-dependent contribution. This procedure gives for the 
rate of energy change 
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The energy current of right movers at the left end of the 
wire is controlled by the noninteracting lead with known 
distribution function P^ . To linear order in V and AT 
a simple calculation gives us 
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It is convenient to combine electric and energy currents 
into the heat current 



JQ= 3E- ^J■J , 



for which we find from Eqs. (HH) and (|27] 
2^2 



3h 



-TAT - 



JQ 



AR 



(29) 



(30) 



where Q^ — E^ — jiN^ is the heat transferred into 
the right-moving subsystem by electron collisions. As 
expected when Q^ = we recover from Eq. (pO| the 
thermal conductance of noninteracting electrons Eq. ([2]) , 
while a nonvanishing Q^ results in an interaction- 
induced change in K. Equation (15(11) was reported earlier 
for uniform wires^^^ and as it is shown here remains 
valid even in the inhomogeneous case. 

So far our basic Eqs. ([22]) and pO|) contain four un- 
known entries: two response currents / and jq due to 
applied bias V and temperature difference AT, and two 
rates N^ and Q^ nonvanishing due to interactions. In 
the next sections we demonstrate that all four quanti- 
ties can be expressed in terms of only two parameters 
u{x) and A/i(x) — fi^{x) — n^ix) which enter the dis- 
tribution function ([9]) of partially equilibrated electrons. 
This would give us the closed set of equations that relate 
{/, Jq} ^ {V, AT} and thus determine the transport co- 
efficients of interest. 



C. Currents I and jq in the partially equilibrated 
wires 

Electric and heat currents can be easily found knowing 
the distribution function f{x,p) of electrons in the wire. 
As we discussed in Sec. |ll] for the partially equilibrated 
wire this distribution is given by Eq. (j9|). It is worth 
emphasizing that for the inhomogeneous case all four pa- 
rameters entering Eq. (|9]): velocity m, chemical potential 
of right- and left-moving electrons ^^^^\ and effective 
temperature T, are, in principle, coordinate dependent. 
Furthermore, the distribution, Eq. (jH]), does not apply to 




FIG. 2: [Color online] (a) Dominant three-particle collision 
which changes the number of right-moving electrons, (b) 
Equilibration mechanism, multi-step backscattering of right 
mover into the left mover. At low temperatures each step 5p 
in momentum space is of the order of '^ T/vf- (c) Energy- 
conserving two-particle scattering process that violates con- 
servation of momentum. This process is possible due to the 
presence of inhomogeneities. 



particles near the bottom of the band, for \p\ < y'mT, 
as explained later in the text (see also corresponding dis- 
cussions in Re f. 1351). This, however, does not cause any 
extra difficulties since transport quantities of interest are 
determined by the behavior of the distribution function 
near the Fermi level. 

By using the definition of the current, Eq. (fT9)) . and 
distribution function (jH]) we obtain the electric current in 
the partially equilibrated wire 



/ 



2e 



Afi{x) + en{x)u{x) 



(31) 



Notice here that although contributions to the current 
due to the electron drift, en{x)u{x)^ and partial equili- 
bration between right and left movers, 2eAfi(x)/h, are 
individually coordinate dependent, their sum must be 
constant along the wire. This is a consequence of the 
particle number conservation. 

Since the heat current jq also does not depend on posi- 
tion, it can be calculated at any point in the wire. In the 
regions not too close to the leads the distribution func- 
tion is expected to have the partially equilibrated form 
(121). Then using expressions (fTO)) and (^5]) for j and je 
we obtain after Sommerfeld expansion of the integrand 
to leading order in T/n <C 1, 



rp2 



3Q = 



6 /i(a;) 



n(x)u{x) 



(32) 



where we introduced ^{x) = ^i—U{x). Since jq is already 
proportional to small u ex y, we were able to replace T 
with T within the linear-response regime. Notice also 
that the particular combination n{x)u{x) / fi{x) defining 
heat current jq must be coordinate independent. 

To make further progress we should elaborate on the 
expressions for the rates N^ and Q^, whose explicit 
forms depend on details of the equilibration mechanism. 
This can be done following the idea suggested in Ref. [Sa 
and we show below how N^ and Q^ can be expressed 
through A/i(a;) and u{x). 



D. Microscopic expressions for A'^^ and Q^ 

First let us identify the leading backscattering mecha- 
nism that contributes to N^ . The most favorable col- 
lisions should involve a maximal number of electronic 
states close to the Fermi points. However, due to the 
conservation of total energy and momentum, collisions 
that change the number of right and left movers can- 
not occur near the Fermi level only, and have to involve 
states deep in the electron band. As was pointed out in 
Ref. .31i the scattering process most important in altering 
the current thus typically scatters two electrons close to 
the Fermi points and one electron at the bottom of the 
band, as schematically depicted in Fig. [5^. It is conve- 
nient to think of this collision as a process in which a 
deep hole, corresponding to the outgoing electron state, 
is backscattered by electron excitations close to the Fermi 
level. These excitations are typically associated with a 
momentum change 6p ~ T/vp <C pf due to Fermi block- 
ing. Let us furthermore characterize this process by in- 
troducing three-particle scattering rate 1/Teee, which can 
be approximated by a constant because the initial and fi- 
nal states both lie at the bottom of the band. Since the 
sign of 6p varies in a random fashion from one collision to 
another the hole performs a Brownian motion in momen- 
tum space. The corresponding diffusion coefficient B can 
be readily estimated. The typical momentum change of a 
hole over time t behaves as (Ap)^ ^ Bt. As we assumed 
the hole changes its momentum by ±T/vf once during 
the time Teee, so we conclude that (Ap)^ ^ {T/vpYt/Teee 
for t >• Teee and thus estimate 



B 



T' 



vIt, 



(33) 



p I eee 



The change AA'^^ in the number of right-moving elec- 
trons over the time t ^ {Ap)'^/B for the segment of the 
wire Aa; is given by the rate t~^ times the number of 
deep holes susceptible to be backscattered. The latter 
can be estimated from the probability to find a left- or 
right-moving hole e~'^ ' "^ and the number of states 
ApAx/h available within the typical momentum range 
Ap ^ y/mT of the backscattering processes shown in 
Fig. [2l3. Taking into account that the scattering of left- 
and right- moving holes both contribute to A^^, but with 
opposite signs, one finally estimates 



AiV« ^ - ^^P^ (e-'^"/^ - e-^'/^ 



A^AxB 

( 

hvrriT^ 



-m/t 



(34) 



with A^ = fi^ — /i-^. A careful calculation based on the 
kinetic equation givesSi 



dx 



2Afj.{x) e 



-tL{x)/T 



h{x) 



4(.) = ^^.(35) 



The expression for the diffusion coefficient B is model 
specific. Our preliminary estimate gives _B ex T'^ in the 
case of Coulomb interaction. 

We continue now with the calculation of the rate Q^j 
which consists of two contributions 



n>R 



Qb + Qp ■ 



(36) 



The first one Q^ is related to the same backscattering 
events that change the number of right-movers A^''^. Ex- 
ploring the fact that both rates are caused by the same 
physical processes it was shown in Refs. [SJ and |3^ that 
there is a relation between A^-'^ and Q^, which we gener- 
alize here for the inhomogeneous case. 



dQ, 



dx 



^2/i(x; 



dN^ 
dx 



(37) 



The logic behind this equation is as follows. The 
backscattering processes transform the unperturbed dis- 
tribution of electrons into the partially equilibrated form 
^ with nonvanishing u{x). The two distributions differ 
most prominently at energies within ~ T of the Fermi 
level. One can thus assume that all the right-moving elec- 
trons contributing to N^ are removed from the vicinity 
of the right Fermi point and placed to the vicinity of the 
left one. Each such transfer reduces the momentum of the 
system by 2pF- The other electrons have to be scattered 
in the vicinities of the two Fermi points to accommodate 
this momentum change. In the special case of three- 
particle collisions, the transfer of electron from the right 
Fermi point to the left one is accomplished in a number 
of small steps with momentum change 6p ~ T jvp , and at 
each step one additional electron is scattered near each 
of the two Fermi points, see Fig. [5]d. As a result of the 
rearrangement of electrons near the two Fermi points, 
the local momentum change Ipp of the backscattered 
electrons is distributed between the remaining right- and 
left-moving electrons, i.e., &p^ -I- &p^ = 2pp. Thus 
the energy of the remaining right movers increases by 
5E^ = vpSp^ whereas that of the left movers decreases, 
SE^ = —vpSp^. Then, the conservation of energy re- 
quires 5p^ = 5p^ = Pp- In the end, the energy balance 
for the right-moving electrons consists of a loss of jjl due 
to the removal of one particle from the Fermi level and 
a gain of 6E^ = vppp = 2fj, due to the redistribution of 
momentum. As a result, for every right-moving electron 
that changes direction, AN-^ = — 1, the right-movers en- 
ergy increases by an amount AE^ = /i, so one concludes 
that E^ = -fiN'^ or equivalently Q^ = -2fiN^. Equa- 
tion (j37p follows from here naturally if one applies the 
same argument but for the segment of wire Ax such that 
rates A^-'^ and Q^ are accounted per unit of length in the 
inhomogeneous wire. 

The other contribution Q^ to the heat transferred by 
right movers in Eq. p6|) is due to scattering processes 
that do not conserve momentum, see Fig. [^t for illustra- 
tion. These two-body collisions are possible only in the 
inhomogeneous case. They do not change the number of 
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right-moving electrons, but do change their energy. It is 
expected that this rate is proportional to the velocity u 
of the electron liquid, Qp ex u. Indeed, two-particle col- 
lisions of Fig. [2t involve a right mover with momentum 
p ~ pf and a left mover with momentum p « —pF- For 
these electrons the drift term pu « ±pFU in the distribu- 
tion function ^ can be absorbed into the temperature 
T, such that right movers can be considered as being at 
an effective temperature T^ w ^^(l -f u/vp) while left 
movers at temperature T^ Ki T{1 — u/vp), to linear or- 
der in uM- According to the general principle of statisti- 
cal mechanics thermalization between these subsystems 
involves the energy flow from "warmer" right movers to 
"colder" left movers that is proportional to the difference 
in temperatures between the two, Qp oc T^ — T^ ex u. 
An explicit microscopic calculation of the rate Q^ done 
in Appendix ([^ gives 



dx 



= —2fi{x 



n{x)u{x) 



(38) 



Here i?in is a scattering length scale associated with these 
momentum-nonconserving collisions. Fig. [2J3, 



T(x) T 
16n{x) ^(x) 



(39) 



where the parameter 



T(x) 



Vn-V. 



2kF{x) 



TThVF{x) 



irhvF {x 

2 

dx 



V, 



2kF{x) 



TrhvF(x) 



(40) 



is expressed through the zero momentum and 2kF Fourier 
components of the electronic interaction potential V. The 
complete rate Q^ is thus given by the sum of Eqs. (|37)) 
and (I55t and we find 



dx 



-2^i{x) 



dN'' 
dx 



2n{x) 



n{x)u{x) 
A„(x) 



(41) 



One should make two important comments regard- 
ing this result. First, the effect of inhomogeneity on 
resistivity of a quantum wire was recently addressed 
in Refs. |32||33| assuming that electrons are fully equili- 
brated and thus described by the distribution function 
([5]). This assumption requires that three-body interac- 
tion processes which change the number of right-moving 
electrons dominate over the momentum-nonconserving 
scattering. However, in a situation where both inter- 
action processes happen on a comparable time scale the 
system is frustrated with a finite value of Afi 7^ and 
thus u ^ Vd — I/en, so that electrons are described by 
the distribution, Eq. ([9|), which we used in our calcula- 
tions. 

To elucidate further the origin of the frustration it is 
important to emphasize that momentum-nonconserving 



scattering resists full equilibration of electrons into a sin- 
gle distribution, Eq. ^. These scattering processes re- 
duce velocity u and thus prevent complete relaxation of 
the difference in chemical potentials A/i. Indeed, since 
the current / is fixed by the external circuit then ac- 
cording to Eq. pip decrease in u implies increase in A/i. 
This effect is opposite to that of the equilibration pro- 
cesses due to three-particle collisions which tend to relax 
Afi. 

Second, the correction to resistance of a nonuniform 
wire was obtained in studies'^^''^'^ by calculating the rate 
of momentum change P^ for right-moving electrons. 
Note, however, that in an inhomogeneous system with- 
out translational invariance momentum is not a good 
quantum number. Indeed, the momentum change of 
right movers due to collisions depends on position and 
thus, is not the same for electrons inside the wire, where 
it was calculated (^i^ than that (actual change) in the 
leads. Our present scheme is free from this difficulty. 
As we show in the next section it is really the rate Q^, 
computed from the well-defined energy exchange, that is 
needed to determine the wire resistance and other trans- 
port coefficients. 



IV. TRANSPORT COEFFICIENTS 

We are set now for the calculation of transport co- 
efficients in an inhomogeneous wire. Indeed, our basic 
Eqs. (1^^ and (pH)) provide electric / and heat jq cur- 
rents as the response to the applied bias V and temper- 
ature difference AT. Interaction effects are captured by 
the rates N^ and Q^ induced by particle collisions that 
encode change in the number of right-moving electrons 
and energy exchange between right- and left-movers, re- 
spectively. These rates are defined by Eqs. (|55|) and (|1T|) 
which still contain the unknown difference between the 
chemical potentials of partially equilibrated right and left 
movers Afi{x) and fiow velocity u{x). Equations (j3ip and 
([5^ are the final ingredients that allow to establish the 
correspondence {V, AT} ^ {I,Jq} and thus find trans- 
port coefficients of interest. 

Technically one proceeds as follows. First, employing 
conservation of currents (recall that I and jq are constant 
along the wire) one can express Afi(x) and u(x) in terms 
of / and JQ from Eqs. (|3ip and p2|) . Second, one brings 
these relations into Eqs. ([55)1 and (PT|) to find the rates 
N^ and Q^ in terms of / and jq, which is possible to 
do in quadratures. Finally, Eqs. ([^^ and ([50)1 define the 
desired correspondence {V, AT} ^ {I,Jq}- Employing 
this procedure we find two independent linear equations, 



2e2 



2n^e 
Sh 



TAT - 



ri) 



ejQ 



6 fi 

12 ri/?-f 7^ 



(42) 
2/ri7I, (43) 



where we introduced the dimensionless parameter, 



ri 



dx 



,-M(x)/T 



(44) 



which quantifies the rate of three-particle processes that 
change number of right-moving electrons [see Eq. ([55]) ] , 
as well as the weighted chemical potentials along the wire, 



dx 
n Jo ^i(a;) 



li{x)e 



-tJ.(x)/T 



— 1 

M^ = — 



dx 



fi^{x)e 



-K^)/T 



ri Jo h{x) 
■y^ = f -^n^{x). 

Jo *-in(2;j 

From Eqs. (|42]) and p3)) we find the resistance 



where 



-^ 



r = ri 



AT=0 



h 
2? 



-^[1 + '^]' 



(45) 
(46) 
(47) 

(48) 



rf^^ 



^T^ + rifi^+j^ 



and Peltier coefficient 



n=^ 



t:^T' 



ri^J. 



AT=0 



6e 



f§T2 + rifi^ + 72 



(50) 



In addition we find the thermal conductance 
K^^ ^-^T3 1 

and thermopower 



18h hIt^ 



TifJ-^ + 7^ 



(51) 



5' = - 



V 

AT 



n'T 



rifi 



1=0 6e 2iT2 + riAi2+72 



Predictably, the Peltier coefficient, Eq. ((5(I)) . and ther- 
mopower, Eq. (|52|) . satisfy the Onsager relation 11 = ST. 
Equations (|48p - ([5^ are the main results of this paper. 
In the following we analyze these general expressions for 
a few modeling examples of inhomogeneities. 



A. Uniform w^ire 

First of all, we perform a consistency check for the case 
of a uniform wire, recently studied in Ref. (35l In the ho- 
mogeneous case U{x) — > all quantities defining R, 11, 
K and S become coordinate independent: ii{x) — >■ £1 
so that JJ? = jf = fj,'^ and n -J> (L/^i)e-^/^. At the 
same time 7 — >■ which is a consequence of momen- 
tum conservation: in a uniform wire two-electron scat- 
tering processes, shown in Fig. [^t, are not allowed. As 



a result, interaction-induced correction to the wire resis- 
tance, Eq. (|49p . reduces to 



tq +ri 



(53) 



where tq = 7r^T^/12/i^. In order to establish a connec- 
tion with the notations of Ref. |35| we invert resistance, 
Eq. (^5]) . with r taken from Eq. ([i^ . to get the conduc- 
tance G = Br^ to leading order in T/n <^ 1 and find 



2e2 



1 



TT^T^ 



L 



12 /i2 L + 4 



(54) 



where following Ref. |35| we have introduced the equili- 
bration length 



•7r2 t2 



(55) 



This result shows that for a long wire L ^ i^q the con- 
ductance saturates to its length independent value, which 
still exhibits noticeable power-law correction in temper- 
ature, 6G = -(2e^/ft-)(7r^T^/12/i^), already mentioned 
in Sec. |lT] [Eq. ([5])]. This saturation of conductance is 
expected, since the electronic system reaches full equilib- 
(49) rium. For short wires, €1 <C i ^ 4q, the interaction- 
induced correction to conductance is exponentially small 
50 = -(2e2//i)(i/^i)e~^/^ and scales linearly with the 
length of the wire. It is worth noting that Eq. ([Ml) is only 
applicable to wires longer than ii . This constraint comes 
from the approximations made when deriving the rate of 
change for the right-moving electrons N^ in Eq. ([55)) [see 
Ref. [3^ for details] . 

The Peltier coefficient and thermopower of a uniform 
wire follow from Eqs. ([501 and (|5^ . and read 



n = ST: 



n^T' 



L 



6e fi L 



(56) 



t-oq 



It shows that II grows from exponentially small values at 
(52) L < £cq to Iloq quoted in Eq. ([8]) at L > 4q. The ther- 



mal conductance behaves very differently though. One 
finds from Eq. (1511) in the uniform limit 



K 



2n^T 4q 
3h L + £c 



(57) 



At L <C 4q one recovers the result, Eq. ^, for nonin- 
ter acting wires, but as the length of the wire grows, K 
is suppressed as 1/L and vanishes for fully equilibrated 
wires, see our earlier discussion presented below Eq. ([7|). 
Equations ([5S|) and ([57]) recover the corresponding results 

of Ref. m. 



B. Two wires in series 

As the simplest prototype of nonuniform system we 
study two uniform wires of lengths Li and L2, with dif- 
ferent densities, connected in series with each other. We 
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ignore here small contribution to transport coefBcients 
coming from 7, which is nonzero only in the near vicinity 
of the junction between the wires and vanishes every- 
where else. Applying then Eqs. P5)) - (|T7| to this setup 
we find for the interaction- induced resistance, Eq. (|49)) . 



TT^r^ 



2/1) 



tlfl 



/^2 eq 



Ml 



M2 



L1L2 



12 



i + i% 



(58) 



be 



K = 



2tt^T 



1 



^h 1 + 1^ + If) 
which is natural generalization of Eq. ([5 



(61) 



C. Wire with long-range disorder 



with chemical potentials /i^ = /i — C/^ and equilibration 
lengths e^c^ = (7r2r7l2Atf)^i*^e^'/^ within each wire i = 
1,2. It is of special interest to consider the limit when 
one wire is infinitely long. Upon taking the limit Li 
the expression for r reduces to 



2rp2 



uJ4 



T2 



T 

Ml 



T 

M2 



/2) 



(59) 



The first term of this formula corresponds to the residual 
resistance of the first wire. It is independent of its length 
Li, which is natural since in the limit Li — >■ 00 elec- 
trons reach full equilibration and the interaction-induced 
resistance should saturate.^'' The other contribution to 

the resistance in Eq. ([5^ is due to the second wire, 

(2) 
which, however, ceases to saturate even when L2 ^ ^oq . 

Thus it remains proportional to the wire length L2 and 
could be much larger than the first term. The absence 
of equilibration in the second wire is rather counterin- 
tuitive. This result forced us to reexamine more care- 
fully continuity equations for electric and heat currents. 
We have found that it is not possible to match both / 
and JQ in the wires while simultaneously imposing van- 
ishing chemical potential difference A/i between right- 
and left-moving electrons. Indeed, let us suppose that in 
the limit Li ^ icq and L2 ^ £eq both wires are fully 
equilibrated, so that A/x = in each wire. The current 
conservation then becomes / = eniVdi = en2Vd2, where 
Vdi,2 is the drift velocity within each wire. According to 
Eq. ([5^ the heat current in this case can be written as 
JQi.2 = -f(7r^/6e)(T2//ii_2) and clearly jqi ^ jq2 since 
Ml 7^ M2- The resolution of this controversy is possible 
only if A/i 7^ at least within one wire even though its 
length exceeds the corresponding equilibration length. In 
Appendix |B] we rederived Eq. ((55)) relying on conserva- 
tion laws only. 

Other transport coefficients do not show dramatic 
changes compared to a single uniform wire and their be- 
havior follows expectedly as a natural generalization of 
Eqs. ([56]) and ([57)1 . We find for the Peltier coefficient of 
two connected uniform wires 



n 



^2r2 



Mi4 



-I- 



M2<: 



6e 1 






L2 

/2) 



(60) 



n saturates to tt'^T'^ /6e^it2) depending on which wire is 
fully equilibrated. The thermal conductance is found to 



We now study more generic models of a nonuniform 
wire. We assume only that disorder variations happen 
on the large spatial scale, kp^ <^ b <^ L, and concentrate 



(X) on the case 7 3> T^. In this case the 7r2T2/12 term 
in Eq. (|49p can be ignored, and the interaction-induced 
resistance of the wire r can be written as 



ri 



riSfi^ 



-7 



ri/i^ -I- 7^ 



(62) 



where we introduced S^-^ 



II- 



/i^. Expression (|62]) is 



applicable to any realization of long-range disorder po- 
tential. We now apply it to two special cases which allow 
simple analytical solution. 



1. Weak disorder 

First, let us assume that amplitude Uq of variations 
in the inhomogeneity potential along the wire is small, 
Uq <C T. It turns out that Eq. ([S^ covers three dis- 
tinct regimes depending on the temperature. At lowest 
temperatures T <^ Ti, where 



Ti 



M 



2 In 



I^F^Fbi^] 



(63) 



three-particle equilibration processes are weak due to ex- 
ponential suppression e~^' "^ of the s catt ering near bot- 
tom of the band.^° In this regime ri5y? <^ riji^ <^ 7^ 
and the resistance of the wire, Eq. (|62p . is given by ri. 
It then follows from Eq. (l44l) that to leading order in 



PiL, pi 



1 



-m/t 



T<Ti, 



(64) 



where pi has the meaning of dimensionless resistivity of 
the wire. It is interesting to compare this result with 
Eq. ()54p obtained for a uniform wire. In the limit L <^ i?cq 
using Eqs. (|^5)) and (1551) we extract from Eq. (l5l)) the 
correction to resistance r ~ {L/(.i)e^'^/^ which coin- 
cides with Eq. ()64p . However, there is an important 
difference in the applicability of this result to uniform 
and disordered wires. In the case of uniform wires r — 
(L/^i)e~^/-^ applies only in the short wire limit L <C ^eq, 
or equivalently at temperatures T -^ T^ ~ (u/ln(L/£i). 
For longer wires r saturates to the length-independent 
value r = n'^T^ /12ii^, and thus always remains smaller 
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than contact resistance, r <^ 1. In contrast, in the case 
of disordered wires the result, Eq. (|M|. does not rely on 
the assumption that L ^ £cq- The crossover temper- 
ature, Eq. (155)) . is controlled by disorder and does not 
depend on L. Thus, although the resistivity pi is small, 
for a sufficiently long wire the total resistance r ^ piL 
can be large, r ^ 1. 

At higher temperatures T ^ Ti (strong equilibration) 
the resistance is given by the sum of two terms. 






II 



(65) 



The contribution of the momentum-nonconserving two- 
body collisions induced by disorder, the second term in 
Eq. (|65l). dominates in the temperature regime Ti ^ 
r<T2, where 



A^ 



In this case the resistance of the wire is given by^^ii^ 



(66) 



P2L, P2 



(T) T 
16n p. 



Ti < T < T2 



(67) 



Here we used Eqs. ([39]) and ([47]), set p^ ^ p^ = /i^ to the 
leading order in Uq <^ T, and (. . .) =_Jq ^(. ■ .) implies 
averaging along the wire. In Refs. |32| and |33| the same 
result for resistivity was derived assuming that electrons 
are fully equilibrated. Here we find that in fact appli- 
cability conditions for p2 are more strict and Eq. (j67p 
dominates only in the temperature range Ti <^T <^T2. 
At higher temperatures the resistance is governed by the 
first term in Eq. (|55]) . 



PzL, pz 



^4t3 {5^P 



18V 



T2 < T < ^ , (68) 



where {Sp^) — (U^) 



since p{x) ^ p — U{x). In 
Eq. ([68)) we introduced thermal conductivity of a wire 
X = (27r2r/3/i)4q- Indeed, according to Eq. ([57]) for 
long wires, L 3> ^cq, we have K = KJh. 

Equation ([68]) shows that resistance r is determined 
by the magnitude of the disorder potential rather than 
its gradients. Similar to the case of two wires in series, 
considered in Sec. lIVBl this feature can be traced back to 
the fact that in the inhomogeneous wire electrons never 
reach full equilibration no matter how long the wire is. 
It is also interesting to note that resistance of the wire 
in the regime of strong equilibration [Eq. ()68p ] can be 
understood from purely hydrodynamic considerations.^^ 



2. Strong disorder 

We now relax the assumption of small variations in 
the amplitude of U{x). Consider the case when gener- 
ally smooth profile of the inhomogeneity potential has a 




FIG. 3: [Color online] Enhanced equilibration due to three- 
particle collisions in the wire segment of length It where 
the inhomogeneity potential is majdmal. A larger value of 
g-t^w/T jjgg^j. X ~ xq favors stronger electron backscattering 



in accordance with Eq. ()35|1 . In contrast, momentum-non- 
conserving two-particle collisions occur throughout the wire. 



well-defined maximum [/q 3> T at some point xq inside 
the wire, see Fig. [3] We assume that such strong fluctu- 
ation of U{x) happens in only one place along the wire. 
This situation is likely to occur in practice due to the 
presence of charged impurities in the substrate and/or 
uneven screening of the nearby gates. 

Since e~^^^^/"^ is largest near xq and dN^/dx ex 
g-M(a:)/T^ see Eq. ([35]) . it is natural to expect that three- 
particle equilibration processes are significantly enhanced 
in the part of the wire where U{x) reaches its max- 
imum. Thus this region of the wire gives dominant 
contribution to resistance at lowest temperatures when 
7^ ^ rip"^ ^ riSp^. Assuming smoothness of li{x) at 
X ~ xo as compared to the sharp e~^*^^^/"^ we can com- 
pute the resistance r ~ rihy applying the saddle point 
approximation to Eq. (|44)) . 



^e-'^"/^ 
h 



It = 



2t:T 
|C/"(xo)| 



(69) 



Here po = p — U{xo), ii — ^i(a;o) and £t is the ther- 
mal length associated with the curvature of the potential 
U{x) near xq. 

In contrast to the equilibration processes, dominated 
by scattering near xq, momentum-non-conserving two- 
particle collisions occur throughout the wire. Resistance 
of the wire is controlled by the latter at intermediate 
temperatures when ri/I^ ^ 7^ ^ riSp"^. For this regime 



TL /T{x)p{x) 



16/io \ n{x) 



(70) 



which is analogous to Eq. ([67]) . The only difference is 
that due to strong variations of U{x) spatial averaging in 
Eq. ([70]) involves not only T(x). 

At higher temperatures, when rip^ 3> riSp^ ^ 7^, 
resistance is again dominated by the scattering processes 
near the top of the inhomogeneity potential. In this 
regime the first term in Eq. (|65p . determined by the 
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amplitude fluctuations of U{x) rather than its gradient, 
gives the leading contribution 



TT^T^ 



24 M^ 4q(a;o) 



(71) 



Here we used Sfi'^/jl'^ — T^/2/iQ found within saddle- 
point approximation from Eqs. (|45l) and (j46p to lead- 
ing order in T/^q <C 1. Equation ([7T|) is analogous to 
Eq. (p5)) . with the thermal length £t effectively playing 
the role of the system size. 



V. SUMMARY 

In this paper we studied the transport properties of 
weakly interacting one-dimensional electrons in the pres- 
ence of inhomogeneities. In this system equilibration 
is strongly restricted by the phase space available for 
electron scattering and conservation laws. The resulting 
equilibration length £eq oc e'^' "^ is exponentially large at 
low temperatures and the partially equilibrated state is 
more likely to be realized than the fully equilibrated one. 
Furthermore, inhomogeneities present in the wire them- 
selves resist equilibration of electrons due to momentum- 
non-conserving two-particle collisions. 

Our main results are expressions (P5 | - ([5^ for the 
resistance, Peltier coefficient, thermal conductance and 
thermopower. We find that the combined effect of in- 
teractions and inhomogeneities can dramatically increase 
the resistance of the wire. For the long enough wire the 
induced correction could be much greater than the con- 
tact resistance of noninteracting electrons h/2e^. This is 
in contrast to the uniform case where interaction-induced 
correction to resistance saturates for L ^ i^q and re- 
mains small as (T/fi)'^ compared to the resistance h/2e^ 
of a non-interacting wire. 

The combined effect of interactions and inhomo- 
geneities is different for thermoelectric coefficients. On 
the one hand, when temperature increases, Peltier coef- 
ficient and thermopower grow from exponentially small 
values, Eq. ©, to H = T5' ~ {tt'^T^ /6e)ijL/'i^). On 
the other, this enhancement is not as dramatic as in the 
case of resistance. Indeed, the difference between the sat- 
urated values of Peltier coefficient and thermopower in 
inhomogeneous wire as compared to that in the uniform 
wires is only in the appearance of the renormalized factor 
71/ /i^ instead of the inverse chemical potential l//i, see 
Eq. ([8]) . Conversely, the thermal conductance of the wire 
decreases due to the equilibration from its noninteracting 
value, Eq. ^, to zero at L >• ^oq- 

The lack of complete electronic equilibration in the in- 
homogeneous quantum wires, which is another central 
observation of our work, warrants additional discussion. 
The notion of equilibration appears naturally since ini- 
tially right- and left-moving electrons entering the wire 
from the left and right lead, respectively, are at different 
equilibria with respect to each other due to the applied 



bias or temperature difference. In the case of weak in- 
teractions three-particle collisions constitute the leading- 
order relaxation process. Although the corresponding re- 
laxation rate is slow (or equivalently relaxation length is 
large) due to the required scattering through the bottom 
of the band, complete equilibration between the right 
and left movers is nevertheless possible in a homogeneous 
wire. Once the length of the wire becomes large such 
that exponential suppression of the equilibration effects 
is compensated by a large system size, L :^ £f.q ^ e^^^ , 
the relaxation of the electron system becomes significant. 
Right and left movers eventually equilibrate to the single 
distribution, Eq. ([5]). When viewed in a reference frame 
moving together with electrons this distribution is simply 
the equilibrium Fermi function. Thus, in the stationary 
frame this is distribution with a boost. 

In the presence of spatial inhomogeneities full equi- 
libration is impeded. This is most transparent 
when momentum-nonconserving two-body collisions are 
present, which unlike three-particle processes favor elec- 
tron distribution function without boost, Eq. Q. As a 
result, the electron system is frustrated due to competi- 
tion between two scattering processes and an intermedi- 
ate distribution, Eq. ([S]), is established. 

Interestingly, lack of full equilibration is a more general 
characteristic of inhomogeneous quantum wires, which 
is a consequence of conservation laws. More precisely, 
in inhomogeneous wires it is generally not possible to 
reconcile full equilibration with conservation of energy. 
Technically speaking, this observation comes from the 
fact that in the fully equilibrated state all currents are 
proportional to the drift velocity. In the linear response 
regime one therefore has only a single parameter Vd to si- 
multaneously satisfy uniformity of particle and heat cur- 
rents along the wire, imposed by conservation laws. In 
inhomogeneous wires this is generally not possible since 
for A/i — the ratio JqJI = 7r^T^/6e/i(x) is not con- 
stant along the wire. Thus the electron liquid must re- 
main in the state of partial equilibration with A/i 7^ 0. 
We have illustrated this point explicitly by considering 
the simplest example of inhomogeneity, namely, a junc- 
tion of two uniform wires with mismatched densities. In 
the general case of a wire with long-range disorder the 
consequence of the partially equilibrated state is that 
interaction-induced correction to resistance of the wire 
is determined by the amplitude of the variations in inho- 
mogeneity potential, rather than its gradients. This cor- 
rection may be large compared to the resistance hjle^ of 
non-interacting wires. 
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Appendix A: Details on the calculation of Qp 

For the inhomogeneous wire we describe the electron- 
electron interaction responsible for two-body scattering 
by its general translationally non-invariant form, 



V{x,x')^V ix-x', 



(Al) 



As a function of its first argument V is assumed to be 
Coulombic in nature and thus short-ranged, with vari- 
ations on the scale of a certain screening length due to 
the nearby gates. The inhomogeneity is captured by the 
second argument with a corresponding variations in V 
on the length scale 6, large compared to both, the Fermi 
wavelength Xp and the range of screening. Our starting 
point for the energy-transfer rate from the right movers 
for the segment of the wire hvp/T <C Ax <C fe is the 
following golden-rule expression. 



AO«(x) 



27r /■ depdtpi de^dtk' 



|V(ep,efc;ep',efc')l^ 



h J (27r)4 

S{ep + efc - ep/ - ek'){ep - ep> + ek' - e^) 
[f^iep){l- f%p,))f^ie,){l- f^{e,,))- 
f%p,)il f%p))f^ie,,)il - f^{e,))] , (A2) 



where the matrix element 



|Vp = |V|| 



|Vi_ 



|V2J 



(A3) 



includes three possible scattering processes of spin-full 
electrons from the initial state (e^, ek) into the final state 
(ep',efe'). These matrix elements will be calculated on 
the basis of semiclassical wave functions, 



1pe,±{^) 



1 



^/hvpiji 



: exp 



4 



dx'^/2r 



Uix)] 



(A4) 

which are eigenstates of the free Hamiltonian 
[—h'^d'^/2rn + U{x)]^pe{x) — £ipe{x) normalized ac- 
cording to J dxtpf^±{x)^p*, ^{x) — 27r(5(e — e') and 
Ve{x) — Pe{x)/m. The subscripts ± refer to the 
right/left branches and we also ignored the backscatter- 
ing wave, since it only leads to an exponentially small 
contribution for kpb ^ 1. Focusing on the states close 
to the Fermi energy, we can simplify the expression for 
the eigenstates, Eq. (jA5|) . of the free Hamiltonian into 



ipe,±{^) ~ ■0M,±(a;)exp 



±i{e-fi) f 
Jo 



dx' 



hvpix') 



, (A5) 



where ipf^^±{x) is obtained from Eq. (IA5|) by setting e = fi. 
This allows us to find, for example, the matrix element 



V|| to first order in the interaction 



V||(ep,efc; £p',ek') 



x+Ax exp 
dX 



0/o"'^-'^^S£i^ 



[hvp{x)]^ 



(A6) 



Here we introduced center of mass X — {x + x')/2 and 
relative y = x — x' coordinates, and kp{X) ~ pp{X)/h. 
Using Eq. ^ and expanding the occupation factors 
f^'/^[e) in Eq. (|A2p to linear order in u{x), and split- 
ting the energy-conserving delta function into two as 
J du!d{ep — Cp' — uj)6{ek — ek' + f^), we can complete e^' 
and ek' integrations and find 

AnR(^\ '^^Tu{x) [ dwdepdek ,^, ,, ,,,2 



4w" 
2^2 



/(6p)(l - f{ep - c.))/(6fe)(l - /(efc + u)) , (A7) 



where /(e) is now the equilibrium Fermi function. The 
corresponding matrix element in these notations reads 



V||(c.) 



^+^%^Vo(X)-V2fe,(X) 



X exp 



-2iuj 



[rivp{x)Y 

^ dx' 

hvp{x') 



(A8) 



The shortened forms Vo and V2kF correspond to the 
zero momentum and 2kp Fourier components of the 
potential V(y,A") with respect to its first variable y 
defined as Vo{X) = JdyV{y,X) and V2fc^(X) = 
/d2/V(2/,X)e^2»fcF(X)a_ ^^ ^-j^jg g^-^^gg ^^ ^^^^ ^^ integra- 
tions can be completed by noticing that 

f{e)-f{e±io) 



/(e)(l~/(e±c.)) = 



oTw/T' 



/ + 00 
de[/(e)-/(e±^)] = ±u;. 



(A9) 



As a result, we obtain following expression for the energy- 
transfer rate: 



Ag5|(x) = 



Tu{x) 
8Trhvp{x) 



x-\-Ax 



Vo{Xi)~V2kAXi) 



TThvp{Xl) 



dXidX2 



VoiX2)-V2kAX2)' 



TThvp{X2) 



-l-oo 



duj 



( uj-^/T V exp (-2*^/xi' Tr^) 



hvp{Xi)hvF{X2) 



(AlO) 



\ sinh -^rp 

Being interested in the temperature range T ^ hvp/b, 
where the exponential is rapidly oscillating, we write 



4w^ 



hvp {Xi ) hvp {X2 
92 



■ exp 



dXidX: 



■ exp 



-2iuj 



~2ia; 



X2 



dx' 



Xi fivp{x') ^ 
^^ dx' \ 



X 



^ hvp{x') 



(All) 
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integrate by parts over Xi and X2 and complete the 
remaining energy integral, which gives 4:Trhvp{x)S{Xi — 
X2). Due to the delta function, one spatial integration 
is thus removed and we find as the final result, Eq. p8l) - 
(|40)) . where all three scattering channels were included. 



the current is determined solely by A/x, while a = 1 
corresponds to full equilibration where / = envd- This 
notation makes it possible to express A/Xi in Eq. (|3T|) in 
terms of a^ as / = 2eA^i/h + ail, which was used in 
Eq. (|B4l) . 

Conservation of the heat current, Eq. 



Appendix B: Series-resistance of two uniform wires 
from the conservation laws 



The resistance of a junction between two uniform wires 
with different densities can be found simply by combining 
conservation laws for the currents with the microscopic 
equations for N^ and Q^. Our starting point is Eq. ([^ 
which is naturally generalized to the case of two wires 
connected in series 



2e'V 



= 1 - eNl^' - eNi 



(Bl) 



where N^'^ correspond to the change in the number of 
right-movers within each segment of the wire. Similarly, 
generalization of Eq. (pO|) for the heat balance reads 



2n^TAT 
3h 



JQ - Q? - Q2 ■ 



(B2) 



Equation (IB1|) defines the resistance, Eq. (j48|). where 
interactions-induced correction is given by 



sTVf eN. 



(B3) 



We use now Eq. psp in the limit of uniform wires to 
express rates 7V/^2 ^^ 



eN, 



irR 



L^ 



/(1-a,)-^, A, = 4''e'''/^ , 1-1,2, 



A, 



by introducing the new quantity Ui — Ui/vd, which has 
meaning of the degree of equilibration. Indeed, a = 
corresponds to the limit of no equilibration, such that 



ejQ 



6 ^ll 



ail 



6 M2 



02/ 



(B5) 



imposes a constraint ai//ii — 012/^^.2 = X^ which must be 
constant. With the help of Eq. (|B4I) the expression for r 
in Eq. (|B3[) can be rewritten as 



hi 

Ai 



L2 
A2 



X Ail 



ii 
Ai 



/^2 



A, 



(B6) 



where the unknown quantity x is yet to be determined. 
The way to find it is from the energy conservation. Re- 
call, that Q^ and N^ are related to each other by 
Eq. ((37|) since they are caused by the same scatter- 
ing mechanism. For the uniform wire Eq. (j37p reads 
Qf = -2fi,N^. As a result, for Eq. ^2^ in the case 
of no temperature bias (AT = 0) we have 

jQ + 2^,lN^ + 2^^2Ni' = 0. (B7) 

Inserting here Eqs. (jB4p and (IBSP one finds 



^T\-2^i,{l-x^^l)^~2f^2{l-X^^2)^=0, (B8) 

D Ai A2 

which allows to find x explicitly. 



2^1 



Ai 



2Ai2 



L2 

\2 



2/x2^ + 2/x2 



(B9) 



A2 



(B4) Finally, using this x ^^ Eq. (jB6p and reexpressing A 



through the equilibration length as Ai = (12^f /7r^T^)£, 



cq 



one recovers Eq. (|58l) . We thus conclude that Eq. ((58|) is 
just a consequence of the conservation laws. 
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